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Anti-ferromagnetic spin chains play an important role in condensed matter and statistical me- 
chanics. Recently XXX spin chain was discussed in relation to information theory. Here we consider 
localizable entanglement. It is how much entanglement can be localized on two spins by performing 
local measurements on other individual spins (in a system of many interacting spins). We consider 
a ground state of antiferromagnetic spin chain. We study localizable entanglement [represented by 
concurrence] between two spins. It is a function of the distance. We start with isotropic spin chain. 
Then we study effects of anisotropy and magnetic field. We conclude that anisotropy increases the 
localizable entanglement. We discovered high sensitivity to a magnetic field in cases of high sym- 
metry. We also evaluated concurrence of these two spins before the measurement to illustrate that 
the measurement raises the concurrence. 

PACS numbers: 03.67.Mn, 03.67.-a, 73.43.Nq, 05.50.-|-q 



Spin chains play an important role in solid state 
physicsiiSiSiii^. For example, inelastic neutron scatter- 
ing on SrCuO^^ can be explained by spin 1 /2 Heisenberg 
chain^. Spin chains are also important for information 
theory. Recently, many-body systems attract lots of at- 
tention in the field of quantum information2*iflii2ii^. Spe- 
cial attention has been paid to the entanglement in these 
systems. Entanglement plays the main role as a physical 
resource for quantum information and quantum compu- 
tation. There is a lot in common between quantum sta- 
tistical mechanics and quantum information theory. The 
role of phase transitions for quantum information was 
emphasized in Ref. il^ Most direct relation between cor- 
relation functions and entanglement was discovered by 
F. Verstraete, M. Popp and J.I. Cirac [we shall use ab- 
breviation VPC] in Ref. % They found that correlation 
function provides a bound for localizable entanglement 
(LE). 

The LE of two spins is defined as the maximal amount 
of entanglement that can be localized on two marked 
spins on average by doing local measurements on the 
rest of the spins [assisting spins]. Here we assume that 
we consider a pure state \<j>) of all these spins. The LE 
has an operational meaning applicable to situations in 
which one would like to concentrate as much entangle- 
ment as possible between two particular particles out of 
multi-particle entangled state. Good examples are quan- 
tum repeater^** and spinotronics^-. Let us consider an 
example of N qubits in GHZ state: 

|Gi/Z>=-i=(|00...0) + |ll...l)) . (1) 

We can measure assisting qubits in |±) basis. This will 
force two marked spins into a Bell state [maximally en- 
tangled state of two qubits]. 

Let us proceed to the formal definition for localizable 
entanglement Eij between two spins marked by i and j. 
Consider a pure state of N spins {(j)) [it is normalized 
= 1]. Every measurement basis specifies an ensem- 
ble of pure states £ = {ps, IV's)}- The index s enumerates 



different measurement outcomes. It runs through 2(^-2) 
values. Here {ips) is a two-spin state after the measure- 
ment and Ps is its probability. The LE is defined as 

=max^p,£;(|V^)) . (2) 

s 

Here E{\ips)) is the entanglement of \ips), characterized 
by concurrence. The concurrence C was suggested by 
W.K. Wootersi^ as a measure of entanglement^. By 
definition it is < C < 1. VPC noticed that it is in 
particular convenient measure for LE. It is important for 
us that concurrence for two qubits state [(/)) = a|00) -I- 
6|01) -|-c|10) coincides with maximum correlation 

C(|(/))) =2|ad-&c|. 

It is difficult to calculate LE. Instead, VPC found 
bounds for LE. The upper bomid comes out of consider- 
ing a global [joint] measurement on all assisting spins. It 
can be related to the entanglement of assistance, which 
is the maximum entanglement over all possible states 
of N spins consistent with the density matrix of two 
marked spins. It was introduced by D.P. DiVincenzo, 
C.A. Fuchs, H.Mabuchi, J.A.Smofin, A.Thapliyal and 
A.Uhlmann"'^^. A simple formula for entanglement of as- 
sistance was found in Ref. Let us denote the density 
matrix of two marked spins by pij . Matrix X is a square 
root of the density matrix: pij = XX'^. The entan- 
glement of assistance measured by concurrence is given 
by the trace norm tr[X'^ {(Ty (g) (Ty)X]. Hence, the upper 
bound of LE is: 




The lower bound on LE is expressed in terms of correla- 
tion functions: 
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The lower bound on LE is based on the following 
observation^: Given a state of N spins with fixed cor- 
relation functions Qq/j between two spins (marked by i 
and j) and directions a and /3, there exist directions in 
which one can measure other spins [assisting spins] , such 
that this correlation does not decrease. Using this obser- 
vation, VPC found a lower bound for LE: 
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(5) 



VPC explicitly evaluated these bounds for the ground 
state of the Ising model and showed that actual value of 
LE is close to the lower bound. 

In this paper we consider the ground state of infinite 
anti-ferromagnetic XXX spin chain at zero temperature. 
We also consider anisotropic version: XXZ chain. We 
calculated the concurrence before the measurement [see 
Appendix A] and compare it to LE. Measurement raises 
the concurrence. 



I. XXX ANTIFERROMAGNETIC SPIN CHAIN 

The Hamiltonian for anti-ferromagnetic XXX spin 
chain can be written as 



Here tri™'' , cry"^' , crl'"' are Pauli matrix, which describe 
spin operators on m-th lattice site. Summation goes 
through the whole infinite lattice. The density of the 
Hamiltonian is a linear function of the swap gate. 

Hans Bethe found eigenfunctions of the Hamiltonian 
of the model in 1931—. The ground state was found 
by Hulten in Ref. 12(1 We shall normalize it to 1. Corre- 
lations are defined as averages with respect to the ground 
state. They are isotropic 
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It took a long time to evaluate correlations functions. 
Nearest neighbor correlation can be extracted from the 
ground state energy^ft. Next to nearest neighbor correla- 
tion was calculated by M. Takahashi in 1977, see Ref. I21I 
Recently it was established that all correlations can be 
expressed as polynomials of In 2 and the values of Rie- 
mann zeta function^ at odd argument a^^i^?i^'^'^^'^i^^i^T . 
These polynomials have only rational coefficients. Third 
neighbor correlation was calculated by K. Sakai, M. Shi- 
roishi, Y. Nishiyama, M. Takahashi , see |2^. These 
results give us the following bounds for localizable en- 
tanglement [LE]: 



Ej^j+i > 0.5908629072 , 

Ej.j+2 > 0.2427190798, 

Ej.j+3 > 0.2009945090 . 

Ejj+i > 0.0346527769. 



(14) 
(15) 
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At large distances correlation functions exhibits critical 
behavior. Asymptotic was obtained inS^iS^: 



There is no magnetization ctq {a ~ x,y or z) 

f7„ = (0|a(^)|0)=O. 

This simplifies the lower bound Q of LE^: 



1 > > 



kl'Vl^'^l 
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Let us now use the explicit expression for correlations 
(dilai^^ ai''M(b) to calculate the lower bound: 



ia«a(^)| 



(-1) 



, (21n|z-j|)V^ 



(18) 



This helps us to estimate localizable entanglement 
asymptotically for two spins, which are far away, i.e. 
|i- j| 00: 
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Even better, but more complicate expression for lower 
bound of LE , can be extracted from the paper—: 
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Here the coupling constant g depends on the distance 
|i — j|. It is defined by: 



7^69 = 2\/2^e'^^+'=|i- 



-J\ 



(21) 



Here = 0.5772 is the Euler's constant and c is a pa- 
rameter [normahzation point]. A good choice for c is 
c = — 1. This boundary for LE is suitable for the full 
range of distances. In Appendix A we calculated con- 
currence before the measurement. It is nonzero only for 
nearest neighbors, see Ref.l38land (|SDJ). Clearly, the mea- 
surement raises the concurrence. 



II. CRITICAL XXZ ANTIFERROMAGNET 
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Let us consider effects of anisotropy of interaction of 
spins. The Hamiltonian of the XXZ spin chain is: 



H' 



XXZ 



E 



,(m)^(m+l) _j_ ^(m)^(m+l) 



+A(a(™Vi"+i) - 1)} . (22) 

We shall consider critical regime (— 1 < A < 1) and we 
use parametrization: 

A = cos(7r?7) , < 7? < 1 . (23) 

Let us remind that the case ry = corresponds to ferro- 
magnetic XXX, which we are not considering here. An- 
other case 7] = 1 corresponds to anti-ferromagnetic XXX, 
see previous section. 

The case ij = 2/3 corresponds to A = —1/2. In this 
case the model admits much simpler solution than generic 
Bethe Anzats, see Ref.'sl'.'ss'.'se . In this case the model 
is super-symmetric, see Ref. ,37^ 

Later we shall see that all these special cases 77 — 
0,2/3,1 have high sensitivity to magnetic field interact- 
ing with spins. 

For general values of rj (I23|l there is no magnetization: 

a„ = (0|a(™)|<^)=O. (24) 
LE is bounded by maximal correlation function: 

(25) 



1>£;,, >max|(0|a«ai^)| 



Correlation functions decays as power laws at large dis- 
tances In |i — j| ^ 1/(2 — 277). Leading terms of correla- 
tions are: 



9 



FIG. 1: F {Eq.^ versus 77 



Many people worked on the subject. Important results 
are obtained inS2. A good collection of references can 
be found in the book [III, see pages 512, 549-553. Since 
< 77 < 1, it becomes clear that ax correlations asymp- 
totically dominate the lower bound: 



UL\ = \Ql'y\ > \Q\ 
Finally we got the following bound for LE: 

E,,>F\i-j\-^ 



(26) 



(27) 



This shows that anisotropy raises the lower bound for IjE. 
The coefficient F was calculated22i2^ 



F = 



2(1-77)2 



V 



2V^r(5^) 



exp 



dt 



sinh(77<) 



t \ sin h(t)cosh((l - 77)i) 



776 



(28) 



The plot is shown on Fig.^ There is singularity in fmic- 
tion F for 77 near 1: 



F ^ (1 — 77) 2 when 77 1 — 



(29) 



Special case 7; = 1 corresponds to XXX antiferromagnet. 

In Appendix A, we evaluated concurrence before the 
measurement. It vanishes at finite distance, see (|86ll . 
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III. XXX ANTIFERROMAGNET IN A 
MAGNETIC FIELD 



Now let us discuss the upper bound (3) of LE. Because 
of translational invariance we have: 



Let us come back to XXX model 
Hlxx = E ^i^^^i"^'^ + + (30) 



(38) 



But now let us add magnetic field 



"-XXX — ^XXX 



(31) 



At large space separations \i 
be simplified 



- i\ cx), correlations can 
T^. This means that both 



s!^ approach (1 — cr^)^- Finally the bounds for LE for 
large 



J I are 



This introduce anisotropy in a different way. In small 
magnetic field ft. — s- 0, small magnetization develops: 
(Jz = h/n'^. For stronger magnetic field, magnetization 
increases. As magnetic field approaches its critical value 
he — 4, magnetization approaches 1 [ferromagnetic state 
with all spins up] : 



(39) 



So magnetic field increases lower bound and deceases up- 
per bound. When magnetic field are close to the critical 
value, the bounds become: 



0-2 = 1 \Jhc-h , 



hr - 0. 



(32) 



A{h)\i 



-]\ 



1 4 / 

" 2 < Eij < —y/hc — h 



(40) 



Exact expression for magnetization az at arbitrary value 
of magnetic field is given on the pages 70-71 of the book 
ITll Averages of other spin components over ground state 
are zero: Ox = Oy = 0. 

In this section, we are considering moderate magnetic 
field < h < he- Asymptotic of correlation functions at 
large distances can be described as follows: 



In Appendix A we evaluated concurrence before the mea- 
surement. It vanishes at finite distance, see (|93|) . 

In the most general case of XXZ a magnetic field cor- 
relations can be described by the similar formulae, but 
parameters (Jz, 9 are different. We shall elaborate in 
the next section. 



1 , . cos{7r{l -az)\i- J \} 

T- ^+"^2 p (OOj 



IV. XXZ ANTIFERROMAGNET IN A 
MAGNETIC FIELD 



Here double bracket in the left hand side means that we 
subtracted ct^: see (4). The coefficients A\ and depend 
on magnetic field. For small magnetic field, critical index 
d is close to 1: 



Let us add interaction with a magnetic field to XXZ 
spin chain: 



H 



XXZ 



XXZ 



(41) 



B ^ l + [2ln{hx/h)Y 



0; (34) 



and for the values of magnetic field close to critical point, 
9 is close to 2: 

1 



Small magnetic field leads to a small magnetization az 
xh. The magnetic susceptibility x is: 



1 



X = 



V 



= 2(1 - -y/hc - h), h^hc] h<hc^A. (35) 

TT 

In Appendix B we discuss the dependence of 9 on mag- 
netic field for intermediate fields. Fig. 4 for 77 = 1 
shows that is a monotonic function of the magnetic 
field. Asymptotic of other correlation functions are: 



7r?7 sm TTTy 



(42) 



A{h)\i~3\ 



-i/fl 



(36) 



Coefficient A{h) vanishes as magnetic field approaches 
the critical value. Exact formula for 9 at any value of 
magnetic field can be found on the pages 73-76 of the 
book 11 and in^^. It shows that 1/2 < 1/6* < 1 < 6^ < 2. 
This means that the lower bound of LE is dominated by 
Gx correlations again 



Here we used parameter 77 related to anisotropy A — 
cos 7777. The dependence of % on 77 is illustrated in the Fig 
2. Let us discuss the plot. The meaning of a singularity 
at 77 = is the following: The case 77 = corresponds to 
ferromagnetic XXX. At zero magnetic field it has spon- 
taneous magnetization pointed in an arbitrary direction. 
Weak magnetic field will align spins to the direction of 
the magnetic field. This makes susceptibility infinite. As 
ry approaches 1, susceptibility approaches I/tt^ [antifer- 
romagnetic XXX case] . 

For stronger magnetic field, magnetization increases. 
As magnetic field approaches its critical value he — 2(1 — 
A), magnetization approaches 1: 



E.j>A{h)\i-j\-^/' . 



(37) 
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0. 



(43) 
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FIG. 2: X (Eq.BU versus tj (Eq.^ 



FIG. 3: ai (Eq.g^, cia (Eq.^ versus r? (Eg.l^ 



Averages of other spin components over ground state are 
zero: ax = (Jy = 0. Here we are considering moderate 
magnetic field < h < he- Asymptotic of correlation 
functions at large distances can be described by the for- 
mulae similar to XXX in a magnetic field case see (|33|l . 

but critical index 6 is different. A formula for 9 
depends on anisotropy A = costtt^. Let us first discuss 
small magnetic field h 0. Critical index is quadratic 
in magnetic field for < < 2/3: 



The coefficient 



ai 



9^- (l + aih^) 
r] ^ ' 



47r?7tan [ 2(1%) ) ^"^^ ^'^ 



(44) 



(45) 



shows singularity at points ry = [ferromagnetic XXX] 
and 77 2/3 [A = -1/2 case]: 



ai 



1 1 



for 77 + 



1 



8l7r2 (?7-2/3) 



for rj 



(46) 
(47) 



The nature of dependence of the coefficient ai is illus- 
trated on Fig 3. In case 2/3 < 77 < 1 small h behavior is 
more complicated: 



(48) 



Notice that the power of magnetic field changes mono- 
tonically from 2atr7 = 2/3to0atry = l. An expression 
for the coefficient a2 is more complicated: 



a2 



2?7e 'I Hq^^ ^ 



tan 



I + i 



r2 



(49) 



Here 



and 



P = (1 — 77) ln(l — 77) + 77 In 77 



4?7\/7r sin 7r77 3 
ho = :—!-e^^^ 



3-2?) 



(^2(^1) shows singularity at point rj ^ 2/3 

for t; - 



(50) 



(51) 



1 



1 



a2 



8l7r2 (7/ -2/3) 



- +0 (52) 



The nature of dependence of a2 on 77 is illustrated on 
Fig. 13 We see that at 7/ = 2/3 critical index 9 strongly 
depends on weak magnetic field. It also depends strongly 
on weak magnetic field for 77 = 1, which is antiferromag- 
netic XXX case. 

For magnetic field close to critical, the index 9 ap- 
proaches 2: 



Trtan [^j tamrri 



(53) 



In Appendix B we discuss the dependence of 9 on mag- 
netic field for intermediate fields. Fig. 4 shows the de- 
pendence of 9 on magnetic field for different values of 77. 
Note that the dependence is monotonic. 

For XXZ in a magnetic field the lower boimd of local- 
izable entanglement is also given by ax correlations 



Eij > 



•lai'^ayh 



(54) 



Asymptotic of the ax correlations is still given by the 
formula (|36|l with 9 described in this section. 



V. SUMMARY 

In this paper we showed that correlations in spin chains 
are important not only for condensed mater physics and 
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statistical mechanics but also for quantum information. 
We considered boundaries for localizable entanglement 
in the ground state of antiferromagnetic spin chains. We 
showed that anisotropy raises the localizable entangle- 
ment. We also calculated concurrence before the mea- 
surement to illustrate that the measurement raises the 
concurrence. 

There are still two open problems left: i. to prove that 
localizable entanglement coincide with the lower bound 
Eij = maxc (IQaal)- ii- to calculate locahzable entan- 
glement for positive temperature. 
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Appendix A 

In this Appendix, we calculate the concurrence be- 
tween two marked spins before measurement. We show 
that the concurrence between i-th and j-th qubits van- 
ishes at finite distance \i — j\ before the measurement. 
The density matrix p of i-th and j-th spins can be rep- 
resented as: 



To calculate concurrence we need 



P = (cTy (Ty)p*{ay Ig) ay). 



(56) 



Here p* is the complex conjugate of p. Subindex /i runs 
though four different values p =^ 0, x, y, z. Pauli matrix 



CTO = / = 



1 
1 



-i 

1 



Following W.K.Wooter^ we define 



1 

1 

1 

-1 



(57) 
(58) 

(59) 



We shall denote the eigenvalues of R by Afe in a decreasing 
order 

Ai > A2 > A3 > A4. 
Then concurrence (C) can be expressed as 
C = max{0, Ai — A2 — A3 — A4}. 



Let us start with the most general case XXZ in a 
magnetic field: 



ai-")^(j and < ct < 1; (60) 



0; (araUi) = (aWai^^)=0;(61) 

{a'^^a^^) = G(* - j) =a^+ gS - j); (63) 
< \9zii ~ j)\ < \9xii ~ j)\ < 1 for large \i-j\.{64) 

Hence density matrix p for i-th and j-th qubits (spins) 
can be expressed as 

9x (^ — J ) 

P = ^ {(7x g CTx + CTy g CTy) 

Gil — ?) (7 

H ^ (7^®(j^ + -{I ®a^+ <7^® I). (65) 

All coefficients ct, gx{i — j), G{i — j) are real, so 

p* ^ p (66) 



P = i<Ty (g) (Jy)p{(Jy ® y) 



(67) 



First three terms in Eq. 1651 commute with Uy ® Uy and 
the last term anti-commute. Let us define 



4 4 
Gi^-J) 



{<Jx ®ax+<7yig ay) 



-a. tit) a. 



and 



m = — (/ (g) CT7 
4 



(55) Notice that [/ ® ct^ -I- CTz (8) /, a^^ax^ ay 



we have 



p = P(i + m, p = pQ~m, [po,m]^Q, 
[p, m] = 0, [p, m] = 0, [p, p] = 0. 



(68) 
(69) 
0. So 

(70) 
(71) 



Now we can simplify the expression for the matrix R: 

R = \NppVp 



'PR 



pI - TV? 



(72) 



Using this representation we can diagonalize R. Corre- 
sponding four eigenvalues {Afc} are: 

V(l + G{^ - j))^ 4ct2, 1 V(1 + G(^ - j))' - Aa\ 

l-G{i- j) gx{i- j) l-G{i- j) _ gx{i-j) 
4 2 ' 4 2 

Now let us consider separately special cases: 

L XXX model with ft, — 0. 
In this case, ct = and 

G{i - j) = gx{i - j) ^ 9z{i- j) = 5, 
The set of eigenvalues of R becomes 

l + .g 1 + g 1+g l-3g' 



73) 



(74) 



To get am explicit expression for concurrence we 
need to consider separately two cases: 
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A: g > (for |i — j| = even) 

1 + .9 ^ l-3ff 



Ai — A2 — A3 



.A4 



4 ' 4 
We can calculate the concurrence C 

C = niax{0,-i-^} = 0, 



(75) 



(76) 



since \g\ < 1. 
B: .g < (for j| = odd) 

Ar = ^^,A.=A3 = A, = i±^ (77) 

The concurrence 

C = max{0, ~^^~^ } (78) 

is non-zero only if g < — i . This happens 
only for j — i ± 1 with 

, , , , 1 - 4 In 2 , , 

5 = = 3 ~ -0.591. (79) 

Hence we have 

1 



Qj+i = ln2 - - w 0.193, 
Cjj+k = if fc > 1. 



(80) 



So, concurrence is non-zero only for nearest 
neighbors [for ground state of XXX model 
without magnetic field]. It was first discovered in 
Ref.S 

II. XXZ model at h = 0. 

In this case (7 = 0. At |« — j| — > 00 

gA^-j)>0, 9xii - j) > \gz{i " j)\ (81) 

and both gx{i — j) and g^^i — j) decay as a 
function of the distance |j — :/|. Eigenvalues of R 
become 



Ai = 



A4 = 



4 2 ' 

X _ X _ i + g4»-j) . 

^^2 — .A3 — , 

i-gz(«-i) 9x{i-j) 



4 2 

Now we can calculate the concurrence: 



(82) 
(83) 
(84) 



C = max{0, ~]^+gx{i~i)~ \9z{i ~ j)} (85) 
It vanishes at distance larger than \i — j\min 

7^ — 9x{\i ~ j\min) — TT Tm 



(2^^)-. 



(86) 



III. XXX in a magnetic field. 



G{i^j) = a^+gS~3) (87) 

For large |« — i| both gz{i — j) and gx(i — j) 
become small. For magnetic field smaller then 
critical 

9z[t~Jj-p — < 1 



(l-(72)2 

\9xii-j)\ > \9z{i-j)\- 
Eigenvalues of R in Eq. [731 become 



(88) 



(89) 



A2-A3-^^ + ^ j-^, (90) 

\ 1 9z{i-i) ,9xii--j). .„,^ 
A4-^ ^ I ^ 1, (91) 

and concurrence becomes 
C = max{0,|5.(^-j)|-^^-^^^i^}. (92) 
Hence concurrence vanishes at distance lager than 

I j\ min 

= \9x{\^-3Un)\^\A{h)\\^~J\-Jl' 



2A(fe) 
l-a2 



(93) 



Appendix B 

In this Appendix, we discuss the dependence of critical 
^ex^ onent 9 on the magnetic field h. We follow the book 

I. Let us start from XXX model. 

Energy of a magnon e(A) is defined by a set of 
equations: 

e(A) - ^ y'^if(A,M)e(M)dM = eo(A), (94) 



-2 



eo(A) = 2/i-^-j-^(95) 



With extra condition e(±A) = 0. This set of 
equation determines the dependence of A on 
magnetic field h. Here A is a value of a spectral 
parameter at the Fermi edge. An important 
object is the fractional charge Z{\): 



(96) 
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The critical exponent is equal to: 
e = 2Z2(A) 



(97) 



For XXX model, the critical field he ~ 4. For 
large magnetic field > he) A = 0. If magnetic 
filed approaches the critical value from below, 
then 



A = i ^Jhc~ \h\ and 



2-1a 



If the magnetic field is small ^ then 



A = — In 
27r 



oo and 6 — \ 



2nA 



(98) 



■ 1 
(99) 



II. Now let us discuss XXZ model. 



In this case we can use the same set of equations 
H94I95|I with K{X,^) and eo(A) replaced by: 



sin(27r77) 



sinh(A — /-t + inrj) sinh(A — /i — iTrrj) ' 



eo(A) = 2h- 



2 sin (Trrj) 



cosh(A+ ^)cosh(A- ^^) 



(100) 



For small magnetic field 



A = (1 - 7]) In 



ho 



oo when /i -> (101) 



Here ho is given by (|51|l . 

But for magnetic field close to critical: 



The critical value of the magnetic field is: 
/ic = 2(1 - A) = 2(1 - cos(7r77)) 



(102) 



(103) 



For general magnetic field h, we solved these 
equations numerically and found that both A{h) 
and 0{h) are monotonic functions of h. The 
numerical solution for 9{h) was shown in Fig 4. 
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